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1. Introduction 

The concept of gluons as elementary excitations in QCD is adequate only at distances 
much less than the correlation length of the theory (the typical hadron radius). Because 
of the dimensional transmutation, the effective gauge coupling constant seizes to be a 
small parameter in the infrared domain and the standard perturbation theory collapses at 
distances comparable with the hadron radius. 

An idea which is perhaps as old as QCD is that infrared- stable elementary excitations 
are described by fluctuating surfaces or strings. The Wilson's strong coupling expansion 
1^ and the expansion of 't Hooft are the two basic arguments in favour of the 
possibility of such a description. 

From the point of view of the expansion, the U{N) gauge theory resembles 
a theory of strings with interaction constant and its multicolor limit — > cxd is 
expected to be described by some kind of a noninteracting string. The formation of strings 
at large distances can be seen as a condensation of dense planar diagrams in the continuum 
theory (fig. 1). However, the only unambiguous evidence for the relevance of strings in the 
perturbation theory is the topology of the Feynman diagrams. Because of the divergencies 
and the need of gauge fixing, it is not easy to make this statement more rigorous. 



Figure 1. A dence planar diagram for the Wilson loop average in the multicolour QCD 

On the other hand, there is the strong coupling expansion in lattice gauge theories 
0]. The gauge theory is discretized so that the gauge field is associated with the links 
£ =< xx' > and the action is a sum over the elementary squares (plaquettes) of the 
lattice. From the point of view of the strong coupling expansion the lattice gauge theory 
is statistical mechanics of random surfaces with contact (self)interactions. In the lowest 
orders of the inverse coupling P = 1/A, the Wilson loop average is given by the sum over 
lattice surfaces spanning the loop. In high orders, the dominant contribution comes from 
overlapping surfaces subjected to strong contact interactions and forming new topologies, 
and the original string picture is lost. 

Some years ago, it has been suggested that the large N limit of the strong coupling 
expansion might correspond to a theory of noninteracting surfaces. In order to study 
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this possibility, the most commonly used character expansion was not the most efficient 
mathematical formalism because it produces diagrams for the partition function and not 
directly for the free energy. In spite of the considerable efforts no evidence for a free 
string picture was obtained 0. 

An alternative formalism was then suggested by V. Kazakov [^. The idea of 0, 
borrowed from the Stanley's solution of the large-A^ vector model, was to reformulate 
the U{N) lattice gauge theory as a kind of Weingarten theory with additional interactions. 
The Weingarten model [§] is a simple (but rather pathological) gauge theory which is 
equivalent to a model of noninteracting random surfaces in the large limit. In the U{N) 
lattice gauge theory one can impose the unitarity of the link variable by means of 
a Lagrange multiplier and replace the Haar measure on the U{N) group by a Gaussian 
measure over complex matrices. The idea of has been accomplished in refs. 0, p!0| ]. 

In ref. |Tl[] we suggested the possibility that the strong coupling lattice surfaces and 
planar Feynman-'t Hooft diagrams appear as two extreme cases of the same perturbative 
expansion. This means that the strong coupling lattice string is perhaps not completely 
irrelevant to the continuum limit of QCD. 



Here we explain the construction of ref. using the standard formalism of the 



mean field analysis in gauge theories [0. Our random surface Ansatz will be defined for 
any choice of the one-plaquette action. In D = 2 dimensions it reproduces the continuum 
limit of the gauge theory in spite of eventual Gross- Witten transition [|13| separating the 
strong and weak coupling domains. However, the first-order phase transition in D = 4 
dimensions is likely to persist and the sum over surfaces has to be modified in the weak 
coupling phase by allowing windows (free boundaries) to appear on the world sheet. In 
the continuum limit the world sheet of the string is tore into a network of thin strips which 
can be interpreted as gluon propagators. 

We start the story with a review of the Stanley's solution of the large-A vector 
model which may help the reader to understand the logic of our construction before dipping 
into the technical details. The analogy with the vector model can help us to get some more 
intuition and avoid the rudest conceptual errors in the case of the gauge theory. 



2. The U{N) vector model 

The fluctuating fleld in this model is an A-component complex vector u^; = 
(u^,...,w^) defined at the points x of the D-dimensional hypercubical lattice and hav- 
ing unit norm 

N 

^ Nevertheless, very interesting results concerning QCD2 have been obtained in a series of 
recent papers |§-|5| using the method of characters. 
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We assume a nearest neighbour interaction /9u* ■ associated with the hnks £ =< xy > 
of the lattice. Then the measure 



AT N 



[du] = Y[ du'du' S{J2 l<P - 1) = du*du S{u* ■ u - 1) (2.2) 



1=1 1=1 



and the interaction are invariant under global U{N) transformations 

^Uu,; lJeU{N) (2.3) 

The partition function 

Z ^ e"^^ = I l[[du^] n e^^"""^ (2-4) 

X £=<xy> 

can be calculated order by order in the inverse coupling constant /3 by expanding the 
exponential as a series of monomials (strong coupling expansion). This leads to a diagram 
technique which in the limit N ^ oo involves tree-like clusters of loops and can be summed 
up explicitly in all orders. 

Alternatively one can transform, by introducing auxiliary fields, the integral over the 
original vector field into a duable Gaussian integral and find an effective field theory in 
which the number plays the role of a Planck constant. 

Introducing a new complex vector field h = (/i*, i = 1,2, ...,A^) we represent the 
unimodular measure as a fiat measure for the unconstrained complex vector field v = 
(1)% z = 1, A^). We insert the exponential representation of the 5- function 

5{u^ -v^) = ^ / rf/i,e^^("'-^') (2.5) 



too 



in the definition ( |2.4| ) of the partition function for every component u^. of the vector field. 
Here hi and are unrestricted complex vectors on which the integration is performed 
with a fiat measure 

N N 

dv = JJ dv^ dv^ , rfh = JJ dh^ dy (2.6) 

Denoting by 



the one-site mean-field integral we write the partition function ( |2.4| ) in the form 
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The potential F{a) is determined by the Ward identity (^^N d / dK^ d / dh* — 
l)e^^ = or, in terms of the a-field, 

r — — + a— - lle^^(") - (2 9) 

which is solved, in the large limit, by 

dF _ Vl + 4a - 1 
da 2a 

(2.10) 



F a = VH- 4a - 1 - In ^- = > /„ — 
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n=l 



Using the integral representation 

giVF(h*.h) ^ J ^^^~^iV[a(h*.h-a)+F(a)] ^2.11) 

we finally write the partition function in a form allowing to perform exactly the integral 
over the vector fields 



Z = e 



NT 



(2.12) 

X l=<_xx'> 



Denoting by C the connectivity matrix of the D-dimensional hypercubical lattice 

D 

(7 = ^^(e^/^"^ +e-^/^"0 (2.13) 

X fl=l 

we can write the result of the gaussian integration over the v, h fields as 

e^^ = / n^"-^"-^'^'^^"^^"''^"^^ ^""^"^"^ = (e'^^°'"^)a (2.14) 

where 

J^Q[a] = -^tr log{l -PaC) (2.15) 

X 

is the sum over vacuum loops of the vector fields. 

Formally the sum over vacuum loops can be obtained by expanding the exponential in 
the last line of ( |2.12| ) as a sum of monomials and applying all possible Wick contractions 

(hrv,,) = (v-h,,) = S,,yS^^ (2.16) 
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The vertex a-ch* ■ h^; is represented by a node at the point x and contributes a factor a^; to 
the weight of the loop. The vertex /3v* Vy is represented by a segment of the loop covering 
the link i =< xy > and contributes a factor /3. The weight of each loop is therefore a 
product of local factors (3 and associated with its nodes and links, correspondingly. 
Denoting by F any closed loop on the lattice, the r.h.s. of ( |2.15| ) can be written as 

^o[/i1=$^/3*''""^n«- (2.17) 
r xer 

The action for the fields ct, ct is proportional to and in the limit N oo these fields 
freeze at their expectation values determined by the classical equations of motion 



or, equivalently. 



dF(a) a/1 + 4q! - 1 

a = ^ = ^?^ (2.19) 

da l-a(3C 

a + aa^ = 1, — ^ = 1 (2.20) 

1 - (3aC 

The last equation has a simple geometrical interpretation in terms of a random walk on 
the lattice. It can be obtained also directly from the unimodularity condition u* ■ = 1. 
Consider the two-point Greens function Gxy = (u* ■ u^). Repeating the same steps as for 
the partition function we find 



^xy 



(v: ■ V,) = (^Ij^b) = J]ci#-^-/3#'-'^^ (2.21) 



(3C 

Therefore the second of the equations ( |2.20| ) is just the condition 



Gxx = 1. (2.22) 

We see that the original vector field u^; acquires a longitudinal degree of freedom and 
instead of 2A^— 1 we have 2N degrees of freedom at each point x. These describe N identical 
noninteracting particles with dynamically generated mass. The constraint u* • = 1 is 
imposed on dynamical level by the condition ( p.22| ) which determines the mass gap of 



the model. After diagonalysing the connectivity matrix we write it as an integral in the 
momentum space 

rr^ = (2.23) 

-™' + EiLi[2sin(p,/2)]2 

= ^-2D (2.24) 

ap 

The continuum limit (infinite correlation length) is achieved when 0, i.e., 2D Pa = 1. 

Let us examine more closely this equation depending on the dimension D. 
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2.1. D <2 

When D <2, the integral is infrared divergent when m ^ and gives 

= / 2I 2 = D < 2 (2.25) 

In the marginal case D = 2 we find the asymptotic freedom law for the dependence 
of the renormalized coupling f3 = /3(m) as a function of the correlation length 1/m 

/3=-^ln^ (2.26) 

2.2. D>2 

Above two dimensions the integral is convergent at m — > and the limit of zero mass 
is achieved at finite value of (3 = (3c. 

In the weak coupling phase (3 > (3c the equations of motion ( p.23|) have no solution 
which means that something in our construction is wrong. This contradiction is resolved 
if we admit that the vector fields acquire nonzero vacuum expectation values (v^;) = Vd, 
(ha:) = hci determined by the saddle-point equations 

hei = /JCvei = 2D/3vei, Vrf = ahrf (2.27) 

which are solved by hd = Vd = ( strong coupling t/(A^)-symmetric phase) or by m? = 
1 — 2D (3a = with hd, Vd 7^ ( weak coupling phase with spontaneously broken U{N) 
symmetry). Geometrically the weak coupling phase is characterized by possibility of the 
the random walk to break into two pieces. Then the two-point Greens function is a sum 
of connected and disconnected parts 

G,y = vl + {x\ \y) (2.28) 

1 — apC 

The unimodularity condition now reads 



Thus in the weak coupling phase there is a spontaneous magnetization Vd = a/1 — Pc/P < 
1 and 2N — 1 Goldstone excitations with zero mass due to the spontaneous breaking of 
the U{N) symmetry. 

The critical value (3c is related to the probability Uq that a Brownian walk returns to 
the origin [O 

2DPc = (2.30) 

i — iln 
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and its large behaviour is given by 

2m=l + ^+3^-| +... (2.31) 

The series diverges aX D = 2 and as L> 2, /3c ^ oo 

2D(3c TFT-^ (2-32) 



TV 



{D-2) 



and below the lower critical dimension D = 2 no transition occurs. In D < 2 there exists 
only a high temperature symmetric phase. 

2.3. D = 1 

The D = 1 vector model has no dynamical degrees of freedom and is analogous to the 
D = 2 gauge theory. 

In terms of the weak coupling parameter A related to by 2//3 = sinh(A/2) we find 

a = tanh — , B = , , , (2.33) 

2' ^ 2sinhA/2 ^ ' 

and the two-point correlator is given by 

G^y = e-(V2)|x-?;| (2.34) 

The two-point correlator can be interpreted either as the sum over all random walks (with 
backtrackings) connecting the points x and y, or as the only minimal path (without back- 
trackings) connecting the two points, with a weight factor Gi = e~^l'^ associated with each 
of its links. In the first interpretation the mean value of a is tuned so that it compensates 
completely the entropy due to backtrackings. In the case of the gauge theory the analog 
of the backtrackings are the folds of the world sheet of the string. We shall see that the 
D = 2 gauge theory can be formulated in terms of minimal nonfolding surfaces. 



3. The U{N) gauge theory 

In the Wilson-Polyakov formulation of the U{N) lattice theory the independent 
fiuctuating variable is a unitary matrix XJi = {{Ui)^, i,j = 1, A^} 

uju, = I (3.1) 
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associated with the oriented hnks I =< xy > of the hypercubical lattice!. The measure 
and the interaction are invariant with respect to local gauge transformations 



U<xy> ^ ^ <xy> 

Uj (3.2) 

The amplitude U(r) associated with the parallel transport along the loop F = {^i^2---^n} 
is given by the product of the link variables 

U(r) = U,,U,,...U,„ (3.3) 

The action in the gauge theory is a functional on the loop fields (|3.3|) . The simplest 
nontrivial loop F is the boundary dp of a plaquette p, an elementary square on the hy- 
percubical lattice. In what follows we will denote by Up the corresponding loop variable, 
the ordered product of the 4 link variables along the boundary of the plaquette p. Let us 
denote by 5'a(U) the one-plaquette action where A is the coupling constant of the lattice 
gauge theory. Then the partition function is defined by the integral over all link variables 



Jees 



where 

N 

[dlJ] = ciUciut5(utu - I) = n dUldiu'f)is[Y.{u'f)^Ui - 6^ (3.5) 

i,j = l k 

is the invariant measure on the group U{N). By JF we denoted the free energy per degree 
of freedom. 

The one-plaquette action is subjected to the following two requirements. First, it 
should be a real function defined on the conjugacy classes of the group U (N). This means 
that it can be expanded in the characters of the irreducible representations of the group 

U{N) 

,iv^s.(u)^ J2 Xn^(I)x,(U)e^^^^^ (3.6) 

— oo<ni <...<njv<oo 

The character of the representation with signature n = {ni, un} depends on the eigen- 
values ui, ...,un of the unitary matrix U as 

X«(U) = "^fJIX^ . (3.7) 



^ We will use the notation i^'^ =< yx > for the link obtained from I =< xy > by reversing 
the orientation. More generally, if Fi and r2 are two oriented hnes, the product rir2 is defined, 
in case it exists, as the line obtained by identifying the end of Fi with the beginning of F2. 
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Second, we assume that the continuum hmit is achieved when A and in this hmit 



S,{V)^-l'^A^ U = e^^^. (3. 



The Fourier image of the action in this hmit is proportional to the second Cazimir operator 

A A ^ 

Sf = ^C,{R) = ^(E-? + - (3-9) 

If we assume that eq. ( |3.9| ) remains true also for finite values of A, the corresponding 
one-plaquette action is known as heat kernel action [|T3[, {U). The exponential of this 
action is a solution of the heat kernel equation on the group manifold 

uX 

The heat kernel action has the nice property of reproducing itself after group integration 

m 

Another simple choice is the Wilson action 5')|^''^(U) = 2^^(U + U^). Its Fourier 
image is given by two different analytic expressions in the limit N ^ oo, depending 
on the value of the coupling constant. For example, for the fundamental representation 
n, = [0,0,..., 0, 1] 



The nonanalyticity at the point A = 1 known as Gross- Witten phase transition is due 
to the fact that the integral for the inverse Fourier transform 



i- I [rfU]e^^^^Wx.(U) (3.13) 



Xn(I) 

is saturated by the vicinity of a saddle point and the saddle-point solution can have two 
analytic branches. In the strong coupling domain A > 1 the eigenvalues of the unitary 
matrix U are distributed all along the unit circle, and in the weak coupling domain A < 1 
the density of the phases of the eigenvalues is supported by an interval [—(j), (j)] with cp < tv. 
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3.1. The gauge theory in terms of nonrestricted complex matrices 

Now let us convert this model into a theory with flat measure and regular interaction 
using the Laplace transform of the measure in the flelds as in the U{N) vector model. This 



has been actually done in the mean fleld analysis of lattice gauge theories [p!2| . We insert 
the exponential representation of the 5-function 

5{ui-vn = ^ dH]e^m-yn (3.14) 

Mil J -ioo 

in the deflnition of the partition function for every matrix element {UiY- of the gauge fields. 
Here and are unrestricted complex matrices on which the integration is performed 
with a flat measure 

N N 

dV = Yl dVldVj, dU= Y[ dWjdH'j (3.15) 

For each oriented link £ =< xy > only one pair of these flelds is introduced and the 
convention 

V<yx> = V^-i = v| = V^^y>, H^y^y = Hg-i = hJ = al^^yy (3.16) 



is assumed. Denoting by 



N'F(U) _ I r ,xTi mr(Htu+Hut) 



[^U]e^^tnw'^+"^'^ (3.17) 



the one-link mean-fleld integral we write the partition function ( p.4|) in the form 

J-JgA^^F(H,)J-J 



^ e 

p 



(3.18) 



where Vp denotes the ordered product of link variables along the boundary dp of the 
elementary square p 

Vp = n (3.19) 

iedp 



The measure is now flat and the integrand is regular. The constraint (|3.1| ) is satisfled in 
the sense of mean values 

((Vjv, -I)...)=0 (3.20) 

In order to give topological meaning of the 1/N expansion we need to express the two 
potentials F and 5"^ in terms of the moments tr( ■ )", n = 1, 2, of their arguments. 
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3.2. The potential for the H-field 

The function F can be determined using the Ward identity 



[tr{d/dU^d/dH) - N]e^^^^^^ = (3.21) 

which is the matrix analogue of (|2.9|) . In the large- limit this function was calculated 
explicitly by Brezin and Gross fl^ in terms of the eigenvalues of the Hermitean matrix 
For sufficiently small field H it can be expanded as a series in the momenta 

a. = ^(HtH)^ (3.22) 

oo 

^H = E E /[.....,M^^^^ (3-23) 

n=l fci ,...,fc„>l 

The coefficients of the series were found in the large limit by O'Brien and Zuber as 
a solution of a system of algebraic relations suggested originally by Kazakov in . These 
relations are equivalent to eq. ( |3.21| ) written in terms of loop variables 



(^^na„_i^+ ^ [(n + /c + l)a„afc + nka^+k-ij^j^] - I) e ^^=0 

n>l k,n>l 



(3.24) 

where we denoted ckq = l,f^n = d/dan- Eq. ( p.24|) is sufficient to determine the 1/A^ 
expansion of the link-vertices /[fci,...,fc„] • In the large limit the identity ( 3.24|) is equivalent 
to a system of recurrence relations 

/[I] = 1 

f[ki,...,k„] + 2^ ^j/[fci+A;j,A;2,...,fc,,...,fc„] + f[kj_-k,L\a]f[k,a] — 

j=2 aePiL) 

where L = /c2, /cn and P{L) is the set of all subsets of L, including the emply set, and 
the bar means omitting the argument below it. The equations for the weights /[^j of the 
disks 

n-l 

/[I] = 1; f[n] + ^ f[k]f[n-k] =0, n = 1, 2, ... (3.26) 

k=l 

are solved by the Catalan numbers 

/n = (-)^"-^^^^^ (3.27) 
^ ' n\{n — 1)1 
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3. 3. The potential for the V field 

The potential S\{y) is defined by the analytic continuation of the r.h.s. of (|3.6| ) from 
the hyperplane V = U, ut U = I to the space of all complex matrices V. 

Using the explicit formula for the characters (|3.7| ) we can express the r.h.s. of ( |3^ ) 
as a series in the moments 

/3. = = ^Vt-; n = 1, 2, ... (3.28) 

The series is convergent for small V and can be represented as exponential of another 
series 

N'SM = N'Y. E (3-29) 

n=l fci,...,fc„/0 



In the continuum limit A ^ the expansion (|3.29|) takes the form i 



oo ^ N _ 



which tends to the 5-function (5(V, I) when V approaches the hypersurface vtv = I. 

Note that the potential ( |3.29| ) for the Wilson action has the form tr(V + vt) only in 
the strong coupling domain. In the weak coupling domain A < 1 it is given by the most 
general expansion ( |3.29| ) with complicated coefficients. Therefore the heat kernel action is 
simpler to work with in the weak coupling regime. 

The lowest-order coefficients for the heat kernel action are, in the large N limit 



S[i]=S[_i]=e / , S[2] =S[_2](1-A)e , 

= _i] = (A - A^/2)e"^, S[i _i] = -e"^ 

3-4- The U{N) gauge theory a la Weingarten 

The interpretation of the functional integral as a sum over surfaces is possible if 
the action is linear in the traces of the matrix fields. This is achieved by the integral 
transformation 

^iv^F(HtH) ^ I f^idaA e^^[%t(HtH)"-H.^]) giv^FM ^3 32) 

n=l 
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This formula resembles eq. (22) in the paper by M. Douglas [19]. We tried to establish an 



exact correspondence between his formalism and ours, but haven't succeed in this 
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By means of another system of parameters we represent the exponential of the plaquette 
action in the form 



n=l 

We therefore introduce at each hnk i and at each elementary square p a set of auxiliary 
loop variables coupled to the moments of the matrix fields 

a^(£),a^(£); n = 1, 2, ... (3.34) 

Pnip)Jnip); n = ±l,±2,... (3.35) 

and represent the integral (|3.18| ) as a theory of Weingarten type described by the partition 
function 



^ (3.36) 

-Q a„(£)tr(HjH,)" "Q [/3„ (p)trV^+/3_„trvi "] \ 

^;n>0 p;n>0 / a,^ 

where the Boltzmann weights of the surfaces are themselves quantum fields and the average 
with respect to them is defined by 

(*)./3=/ n c^a.(£)da.(£)e-^^^^^^+^K^)l 

_ (3.37) 

n rf^.(p)ci/3.(p)e-^^^^^^+^[^(-)l* 



The integration over the matrix fields in ( J3.36D will produce an effective action which 
will be interpreted in the next section as a sum over closed connected surfaces on the 
lattice. The leading term in this action is proportional to A^^. Therefore, in the limit 



N ^ oo the integral over the scalar fields ( |3.34| ) and ( p.35|) is saturated by a saddle point 



and these fields freeze at their vacuum expectation values 

(a^(£)) = 0M) = 0-n{p)) = /3n (3.38) 

This is natural to expect since the auxiliary scalar fields are coupled to gauge invariant 
collective excitations of the original matrix fields (the moments i-)^) which, due to selfav- 
eraging, become classical in the large limit. The classical values of the fields an can be 
determined from the stationarity condition for the effective action, or, alternatively, from 



13 



the Ward identity ( 3.2C1| ). Thus in the hmit ^ oo we obtain a theory of noninteract- 



ing planar {G=0) random surfaces with Boltzmann weights determined dynamicaUy. This 
picture is conceptuaUy the same as in the U{N) vector model. 

Since we are interested only in the large limit, the fields ( p.34|) and (|3.35|) can be 



replaced with their vacuum expectation values even before the integration over the matrix 
fields. We therefore find the following expression for the free energy of the U{oo) gauge 
theory 

(3.39) 

-Q a„tr(HtH,)" "Q g^/3„ (trV^+trv|") 
£-n>0 p;n>0 

The values should be kept as free parameters. The parameters are determined 

by the condition of unitarity ( |3.2C1| ) and (3^ can be considered as coupling constants defining 
the plaquette action. 



4. Branched random surfaces 

4.1. Feynman rules in the large N limit 

Assuming that the saddle point for the integral ( |3.18| ) is the trivial field H = V = 0, 
we can express the free energy as a sum of vacuum Feynman diagrams. 

The Feynman rules are obtained as follows. Inverting the bilinear part of the action 
we find the propagators 

{iHi)){v})f) = {{Hj)]{V,)f) = ^616^; {{v})]{V,)f) = ^616^ (4.1) 

The vertex /3^trVp is associated with the loop [dp)"^ going n times around the boundary dp 
of the elementary plaquette p. We attache to the loop {dp)^ a little planar surface Pn called 
n-plaquette. The n-plaquette pn has area n (we assume that the elementary plaquette has 
unit area) and a boundary going n times along the boundary of the plaquette p. The first 
three multi-plaquettes (n = 1,2,3) are shown in fig. 2 where for eye's convenience we have 
displaced the edges of the boundary. 

The n-plaquette with n > 1 has a singularity (branch point of order n) containing 
Gaussian curvature 2Tv{n — 1) at some point, say, at its centre. The n-plaquettes will serve 
as building blocks for constructing the world sheet of the lattice string. The Boltzmann 
weight of an n-plaquette in /3„. The propagator V — V has simple geometrical meaning: 
it glues two edges of multiplaquettes having opposite orientations. 
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The absence of a. H — H propagator means that the H field just play the role of glue 
for identifying the edges of the multiplaquettes. Besides the simple contraction V — V 
there are cyclic contractions of n pairs of oppositely oriented edges, n = 2,3,..., which 
are represented by the vertices The edges are glued half-by-half in a cyclic 

way. The surface obtained in this way has a singularity at the middle of the link where a 
curvature 27r(ri — 1) is concentrated. We will call as before this singularity a branch point 
of order n. The Boltzmann factor associated with a cyclic contraction of order n is an{i)- 
The simple contraction and the first nontrivial ones are shown in fig. 3. 



n = l n = 2 n — 3 

Figure 2. Multiplaquettes of orders 1, 2, 3 



n= 1 n = 2 n = 3 

Figure 3. Cyclic contractions of n edges, n=l,2,3 

The vacuum Feynman diagrams are surfaces composed of multiplaquettes glued to- 
gether along their edges by means of cyclic contractions. An important feature of these 
surfaces is that they can have branch points (that is, local singularities of the curvature) 
not only at the sites, as it is the case in the original Weingarten model, but also at the 
links and plaquettes of the surface. These singularities can be interpreted as processes of 
splitting and joining of strings. The weight of a closed surface is a product of the weights 
of its elements times a power of N which, with our normalizations, is equal to the Euler 
characteristics of the surface. 
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The free energy of the U{oo) gauge theory is given by the sum over aU closed connected 
surfaces with the topology of a sphere. 

The connection with the traditional Feynman-like diagrams is by duality. Sometimes 
it is mire convenient to use the traditional diagrammatical notations. Then a surface 
bounded by a loop F = [£1^2 •••^n] will be represented by a planar Feynman graph with n 
external legs dual to the links £1, The cyclic contraction of 2n is represented by a 

vertex with 2n lines as is shown in fig. 4. 



Figure 4. Vertices dual to cyclic contractions 

Applying these Feynman rules to the Wilson loop average VF(r) = (|^V(r)) we find a 
representation of VK(r) as the sum over all planar surfaces bounded by the contour F. The 
Boltzmann weight of a surface is the product of the mean values (|3.38 ) associated with the 
multiplaquettes and cyclic contractions. The local Boltzman weights can be decomposed 
into factors contributing to the total area of the surface and the length of its boundary, 
and factors related to the branch points. We denote 

= = ci>(i), (4.2) 

Then the string path integral for the Wilson loop reads 

Vr(F) = e-"^o^P! J2 e-^o^[-51n[cS] (4.3) 
S:ds=r 

where L[F] is the length of the contour F, A[S] is the area of the surface S, the bare string 
tension Mq and the boundary mass mo are related to the original variables cti, Pi by 

-Mn -2 3/ 



e-^^^» = ai/3[ e-"^° = V^i (4.4) 

and the 0-factor is the product of the local factors associated with the branch points. 

The Q-factor of the surface S containing n-plaquettes, n = 1,2,..., and A'l^'* 

cyclic contractions of order n , n = 1,2,..., is given by 



m = n n (-^i^^)^"^' (4-5) 



n=l (j=l,2 

ant the area of the surface is 

A^Yl ^^-'^ (4.6) 

n=l 
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4-2. Irreducible surfaces (a miracle) 

Unlike the random walk, the two-dimensional surface can exist in configurations with 
very uneven intrinsic geometry. The typical syngularity is a "neck" representing an inter- 
mediate closed string state of small connecting a "baby universe" with the main body of 
the surface. The most singular configurations with necks can be readily removed from the 
path integral of the string. Their contribution can be taken into account by modifying 
the local Boltzmann weights associated with cyclic contractions of edges. Miraculously, 
the new Boltzmann weights become easily calculable and do not depend neither on the 
dimension D nor on the choice of the one plaquette action! This is why we will restrict the 
sum over surfaces to irreducible ones which will render the random surface Ansatz much 
simpler than it its original version. 

The configurations to be excluded are surfaces with necks occupying a single link (fig. 
5). We call such a surface reducible with respect to this link.. 



Figure 5. A reducible surface 

In what follows by sum over surfaces we will understand a sum over irreducible surfaces 
defined as follows. 

Definition: 

A surface S is reducible with respect to given link £ e £ if it splits 
into two or more disconnected pieces after being cut along this link. A 
surface which is not reducible with respect to any £ e £ is is called 
irreducible. 

4-3. Evaluation of the weights of the cyclic contractions 

Consider the Wilson loop average for the contour F = (£i^2---) with coinciding end- 
points. 
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We will exploit the unitarity condition vvt = I, applied to the loop amplitude. It 
means that the Wilson average VF(r) will not change if a backtracking piece ££~^ is added 
to the contour F 

W{TU-^) = W{T) (4.7) 
The sum over surfaces spanning the loop rii~^ can be divided into two pieces 

W{TU-^) = W{T)W{U-^) + Wconn{TU~^) (4.8) 

The first term is the sum over all surfaces made of two disconnected parts spanning the 
loops r and The second term contains the rest. The constraint ( [4.7| ) is satisfied for 

all loops if W{£i~^) = 1 and Wconn(r^^~^) = 0. But there is only one irreducible surface 
spanning the loop £i~^] it contains a single contraction ai{i) between i and £~^. Therefore 
CKl = 1. 

Now consider a surface contributing to the second term Wi. The links i and may 
be connected to the rest of the surface by means of the same contraction or by two cyclic 
contractions The condition that their total contribution is zero is 

n-l 

CKn + ^ a/cttn-fc = 0, n = 2, 3, ... (4.9) 

fc=l 

This equation is illustrated in fig. 6 by means of the standard graphical notations (fig. 4). 
Eq. ( ^^ ) is identical to the loop equation satisfied by the coefficients fn in the expansion 
( ^.23| ), which is solved by the Catalan numbers 

«n = /N = H^'^-^^^^^ (4.10) 
^ ' n! n — 1)! 



Figure. 6. Graphical representation of the unitarity condition 
The second derivation is based only on the the sum over surfaces for the trivial Wilson 
loops 

= (V^vJ)-) = 1, n = l,2,... (4.11) 

For each of these Wilson loops the sum over irreducible surfaces contains only finite number 
of terms, namely, the link- vertices contracting directly the edges of the loop £i~^. For 
example, 

W{ir^) = ai, W[{£r^f] = 2a\ + aa, ... (4.12) 
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Introducing the generating functions 

oo ^ oo 

«;(t) = 5^rw^[(£ri)"] = -— , /(t) = i + 5^r(aO'^ (4.13) 



i-t 

n=0 n=l 



we easily find the relation |11 



w{t) = f[tw^{t)] (4.14) 
which is solved by the function generating the Catalan numbers 



. i±4±l ,4.15) 

In this way, unlike the U{N) vector model, the classical values of the auxiliary fields 
do not depend on the dimension of the space-time. The expectation value a in the vector 
model is determined by the long wave excitations of the random walk. Here in the gauge 
theory, the expectation values of the fields an are determined in purely local way, as it is 
clear from their derivation. One way to explain this difference is the local character of the 
U{N) invariance in the case of the gauge theory. 

On the contrary, the weights of the multiplaquettes will depend on the dimension as 
well as on the choice of the one-plaquette action. If wee are interested only in the large 
N limit, it is more convenient not to try to calculate them by solving the equations of 
motion but just to take them as independent coupling constants. The continuum limit (if 
it exists) then will be achieved along a trajectory 

Pn = Pn{X,D), A^O (4.16) 

where A is the coupling constant. 

4-4- String representation of the Wilson loop average. Resume 

Let us summarize what we have achieved by now. The Wilson loop average VF(r) 
in the [/(oo) gauge theory defined on a D-dimensional lattice is equal to the sum of all 
planar irreducible surfaces having as a boundary the loop F. These surfaces are allowed to 
have branch points of all orders located at the sites, links and plaquettes of the space-time 

(k) 

lattice. Introducing unified notations oon for the weights of the branch points associated 
with the /c-cells of the space-time lattice (sites are 0-cells, links ate 1-cells, and plaquettes 
are 2-cells) where 



UJ 



(0) = 1 

n 



^n' = /[ 



{2n - 2) 



n!(n-l)! (4.17) 



13. 



we can write the sum over surfaces as 

W{T) = J2 ^{S)e-^°^''''^^^ (4.18) 
ds=r 

where the factor 0(5) is a product of the weights of all branched points of the world sheet 
and Mq = — In /?i . 
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4-5. The contribution of the surfaces with folds is zero 

The weights of the branch points provide a mechanism of suppressing the backtracking 
motion of the strings or, which is the same, world surfaces having folds. The surfaces with 
folds have both positive and negative weights and their total contribution to the string 
path integral is zero. We have checked that in many particular cases but the general proof 
is missing. 

It is perhaps instructive to give one example. Consider the surface in fig. 7. 



Figure 7. A piece of surface having a fold 

It covers three times the interiour of the nonselfintersecting loop C and once - the 
rest of the lattice. There are two specific points on the loop C at which the curvature has 
a conical singularity; they can be thought of as the points where the fold is created and 
annihilated. Each of these points can occur either at a site or at a link (in the last case it 
is associated with a weight factor uj2^^ = /p] = — 1). 

Let us evaluate the total contribution of all irreducible surfaces distinguished by the 
positions of the two singular points. Denoting by no and ni{— no) the number of sites 
and links along the loop C, and remembering that the a singular point located at a link 
has to be taken with a weight f[2] — —1, we find that the contribution of these surfaces is 
proportional to 

[^^^^^ + /p,«on. + /gi^^i^^] - Im + 2/p,„:l = (4.19) 

The second term on the left hand side contains the contribution of the reducible surfaces 
which had to be subtracted. A reducible surface arises when the two points are located 
at the extremities of the same link (ni configurations) or when one of the points is a 
branch point and the other is located at one of the extremities of the same link (2ni 
configurations) . 
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4-6. Loop equations 

In ref. we proved that the loop equations in the Wilson lattice gauge theory [^| 



are satisfied by the the sum over surfaces ( |4.18| ) . Let us only write here the general formula 



which is derived in the same fashion. For any link £ G F 

oo 

E E Pn[{wir{dpr)-w{r{dp)-^)] = Y,wira')wir,,,Mej')-6ie-\i^^^ 

n=lp:dp3(. (.'eV 

(4.20) 

The sum on the l.h.s. goes over the 2[D — 1) plaquettes adjacent with the link I and the 
closed loops Vii' , F^/^ in the contact term are obtained by cutting the links £ and I' and 
reconnecting them in the other possible way. 

^.7. The trivial D = 2 gauge theory as a nontrivial model of random surfaces 



The case D = 2 was considered recently in details in |21]. For the heat kernel action 



(2) 

the string tension Mq and the weights uJn are given by 



Mo^^; .f^i:(„:O^MrHi-^A..., (4.2, 

m=l 



Let us illustrate how the string Ansatz ( [4.18| ) works for the simplest nontrivial example 
of a Wilson loop (fig. 8) which have been calculated previously using the Migdal-Makeenko 
loop equations p2|. 



Figure 8. A contour with the form of a flower and a branched surface bounded by it 

The nonfolding branched surfaces spanning the loop will cover the three petals of the 
flower (denoted by 1 in fig. 8 ) once, and the head (denoted by 2) - twice. Therefore the 
total area will be always A = Ai -\- 2A2, where Ai is the area enclosed by the domain 
i {i = 1,2). Each of these surfaces will have a branch point located at some site, link or 
cell of the overlapping area 2. 
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Let us denote by no,ni, n2 the numbers of sites, links, cells belonging to the domain 
2. Then the sum over the Q-f actors due to the branch points reads 

J2^= Yl nkuji'"'^ = {no - ni + {1 - X)n2) = 1 - n2X (4.22) 

fc=0,l,2 

(Here we used the Euler relation no — ni + n2 = 1-) Therefore 

W{C) = (1 - ^2)6-^(^1+2^2) (4 23) 



5. Strong versus weak coupling 

Let us first discuss the case of the Wilson action where the mean field problem is 
solved exactly in the large N limit |T^]. It is well known that when > 4 and D = 4, 



the strong coupling phase of a theory with Wilson action is separated from the continuum 
limit by a nondeconfining first-order phase transition. The assumption that Hd = Vd = 
is a local minimum of the free energy is justified in the strong coupling phase of the gauge 
theory. In the weak coupling phase the matrix fields will develop vacuum expectation 
values. The classical fields form an orbit of the gauge group 

V<,j,> = U,VciU-\ U^^yy = U.HciU-^ (5.1) 

The diagram technique of sect. 2 then has to be modified according to the general rules 
explained in fl^ . The surface elements will contain free edges, and, as a consequence, win- 



dows will appear spontaneously on the world sheet of the string. Whether these windows 
destroy or not the world sheet in the continuum limit is a dynamical question. Believing in 
the confinement, we expect that in D < 4 dimensions the windows are still not sufficiently 
large for that, and the world sheet will have in the continuum limit the structure of a dense 
network of thin strips separating the windows (fig. 1). These strips will correspond to the 
gluon propagators in the standard Feynman rules. At the critical dimension D = 4 the 
effective string tension of the string (with windows on the world sheet) should scale with 
the coupling A according to the asymptotic freedom law, just as the mass in the vector 
model does in D = 2 dimensions. 

In D > 4 dimensions the world sheet of the string is eaten by one or several large 
windows and the Wilson loop behaves as exponential of the length of the contour. 

Whether this picture is true in geleral, can be decided by studying the mean field 
problem for a general one-plaquette action. 
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